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Abstract The Hopf algebra of undecorated rooted trees has tamed the combinatorics 
of perturbative contributions, to anomalous dimensions in Yukawa theory and scalar <p z 
theory, from all nestings and chainings of a primitive self-energy sub divergence. Here 
we formulate the nonperturbative problems which these resummations approximate. For 
Yukawa theory, at spacetime dimension d — 4, we obtain an integrodifferential Dyson- 
Schwinger equation and solve it parametrically in terms of the complementary error func- 
tion. For the scalar theory, at d = 6, the nonperturbative problem is more severe; we 
transform it to a nonlinear fourth-order differential equation. After intensive use of sym- 
bolic computation we find an algorithm that extends both perturbation series to 500 
loops in 7 minutes. Finally, we establish the propagator-coupling duality underlying 
these achievements making use of the Hopf structure of Feynman diagrams. 
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1 Introduction 



In JT|, H we developed the Hopf-algebra techniques of || |, |], |] to tame the combinatoric 
explosion of perturbation theory. For Yukawa theory at spacetime dimension d = 4, and 
also for a scalar 3 theory at its critical dimension d = 6, we resummed all nestings and 
chainings of a divergent self-energy skeleton, using the Hopf algebra of undecorated rooted 
trees, to progress beyond the rainbow |7|, || and chain || |l(J approximations for anomalous 
dimensions. Pade-Borel resummation to 30 loops gave stable results at a Yukawa coupling 
g = 30, enabling comparison with resummation methods. However, our previous work 
left unanswered two pertinent questions. 

1. What are the nonperturbative problems whose perturbations [|l], were developed? 

2. Are there limits on the coupling strength for which resummation may be performed? 



Here we answer both questions. Our results provide a stringent testing ground [[□]] for con- 
fronting resummations of symbolically computed perturbation theory with high-precision 
numerical analyses of the corresponding nonperturbative integrodifferential equations. 

Sect. 2 derives the nonperturbative problems as integrodifferential Dyson-Schwinger 
equations. These may be transformed to nonlinear differential equations, of second order 
in the Yukawa case and fourth order in the scalar case. The anomalous dimension is then 
defined, nonperturbatively, by the unique value of the derivative of the renormalized self- 
energy at the subtraction point q 2 = fi 2 that guarantees a well-defined self-energy at all 
euclidean momenta with < q 2 < fi 2 . For the Yukawa case a parametric solution is ob- 
tained in terms of the complementary error function, erfc, with the anomalous dimension 
7 at coupling a = (g/4ir) 2 given by the implicit condition 

= ex P<>o) er fc(po) ; Po ■= ^7= (1) 
V 2tt V2a 



which we solve by a Newton-Raphson method. The asymptotic perturbation series 



( 2 ) 

n>0 Z 

has a coefficient at n + 1 loops given recursively by 

n 

Gn+l — finfi + n 2-s GkGn+l-k ■ (3) 
fc=l 

The corresponding problem in the scalar theory entails boundary conditions on a nonlinear 
fourth-order differential equation. We show how to develop the perturbation series 

n>0 U 

for the scalar anomalous dimension 7, far beyond the 30 loops obtained in After 
attempting a variety of methods, we found one that delivers 500 loops in 7 minutes, 
giving the 1675-digit integer G 50 o of Table 1. 
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In the Yukawa case, an obstacle exists at the critical coupling constant g cr i t : = 2(2n) 3 ^ 2 m 
31.5, corresponding to a critical anomalous dimension 7 crit := —2. The asymptotic series 

is obtained as a — > 0. It is clearly vacuous as 7 approaches —2. At such strong couplings, 
one should use our new nonperturbative result flU), with 

erfc po -=^= d P exp(-p = 1 - — ^ —— — — . 6 

v^r Jp v/tt ^ n! (2n + 1) 

Having thus solved the Yukawa problem at all couplings, we use its perturbative for- 
mulation to motivate methods for the more severe fourth-order scalar problem, whose 
Dyson- Schwinger equation indicates an obstacle at 7 = — 1. 



2 Dyson-Schwinger analyses 

In Yukawa theory, with an interaction term gipuip at d = 4, we consider the dimensionless 
renormalized self-energy term S(g 2 ) in the inverse propagator ^(1 — £(g 2 )) of a massless 
fermion field ip at euclidean 4-momentum q. The subtraction is performed at q 2 = fi 2 , 
with S(/i 2 ) = 0. Then 7 is the value of dlog(l — £(g 2 ))/<ilog(g 2 ) at q 2 = /i 2 . In the scalar 
theory at d = 6 we consider an interaction gftarf) of the massless neutral scalar field a 
with a massless charged scalar field 0, whose inverse propagator is g 2 (l — S(g 2 )). 

The infinite set of Feynman diagrams which we resum comprises all those whose sub- 
divergences result from nestings and chainings of the one-loop self-energy skeleton. A 
12-loop example is given in Fig. 1, with a subdivergence structure encoded by the undec- 
orated rooted tree of Fig. 2. The Hopf algebraic methods of @] gave the integer coeffi- 
cients G n and G n of the perturbation series (QH) up to n = 12 loops, with a coupling 
a := g 2 /(A-K) d l 2 . In |l| we were able to extend this analysis to 30 loops, obtaining 

G 30 = 10272611586206353744425870217572111879288 (7) 
G 3 o = 19876558632009586773182109989526780486481329823560105761256963720 (8) 



each of which subsumes 4.6 x 10 20 subtractions in the BPHZ formalism flT2" |. At that point, 
our dimensionally regularized method hit a limited imposed by memory requirements, 
since we had to store a triple series in powers of the coupling a, the logarithm log(g 2 /yU 2 ), 
and a dimensional regularization parameter e. 

We shall greatly improve on the previous analyses Jl], [|, by using Dyson-Schwinger 
methods. First we show how to reduce the triple series to a double series, by dispensing 
with dimensional regularization. Then we describe the symbolic computation proof meth- 
ods that led to discovery of a single-series method, which dispenses with the logarithm. 
They culminate in a propagator-coupling duality which we then prove thanks to the Hopf 
algebra structure. 
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Figure 1: A twelve-loop example of the type of diagrams considered. 




Figure 2: The corresponding tree structure which determines the necessary subtractions 
to get local counterterms. 



2.1 Integrodifferential equations 

Modulo divergences, the Yukawa problem corresponds to the Dyson- Schwinger equation 
g 2 X(g 2 ) = / - — f„~ ,, ON , q L - subtractions (9) 



7T 



with subtractions that give S(a* 2 ) = 0. 

Let us pause for a moment and consider the structure of this equation. Its design 
guarantees that it represents a sum over all undecorated rooted trees: If we let 



q 2 ll{q 2 )[t] = ^ j dHV{l 2 )[B_{t)] 



i 2 (g + iy 



be, for each q 2 , a character on the Hopf algebra of undecorated rooted trees (which implies 
to define E(<? 2 )[e] := 1 and E(g 2 )[tit 2 ] := S(g 2 )[ti]S(g 2 )[t 2 ]), then we can regard 

D + — I <Vl- q ' 1 



n 2 J l 2 (q + l) 2 

as the corresponding representation of the B + operator with coproduct A(T) = T £g> 1 + 
(id®£+)A(B_(T)) §, 



5 + (S)(g 2 )[t] := nq 2 )[B + (t)] = ^J d^(l 2 )[t]- 



q-l 



'/ 2 (g + /) 2 ' 

The Dyson-Schwinger equation can then be written as 

q-l 



g 2 S(g 2 )[X] = ^/rf 4 /S(/ 2 )[l/(l-X)] 



/ 2 (g + /) 2 ^ 
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where the series in Hopf algebra elements X is defined by X = aB + (l/(l — X)), a form 
which we will use later on. 

As we are confronting a character of the Hopf algebra of undecorated rooted trees, 
prominent in the study of differential equations in the context of Runge-Kutta methods, 
we expect to be able to find a solution to the Dyson-Schwinger equation turning it into 
a differential equation. This is indeed the case. We proceed by first turning the problem 
into a integrodifferential equation. 

The integration over the 4-dimensional euclidean loop momentum, /, factors into an 
angular and a radial part. The angular integration of 

,■1 1 , 2 + / 2 (1Q) 



(q + l) 2 2 2(q + l) 2 
is performed by using the 4-dimensional angular average 



WtT 1 ^) =min(g 2 ,/ 2 ). (11) 

d=A 



We split the radial integral over y := I 2 into parts with y < q 2 and y > q 2 . For the latter we 
fix the subtractions by requiring that S(/x 2 ) = 0, obtaining the Dyson-Schwinger equation 



= \ %7-r + iV) " *V) (12) 

1 V y - yZ{y) 

- Sr^7Tf?V- (13) 



2 Jo y-yY,{y) \x 
Differentiating w.r.t. x := q 2 we obtain the integrodifferential equation 



3 d w_a r y d y 



xVS(x) = a / — ¥—Jr—— (14) 
dx Jo l-E(y) 

which makes no reference to the subtraction point. 

To repeat this analysis in the scalar case, we use the 6-dimensional angular average 

which gives the Dyson-Schwinger equation 

r« 2 ciy 



^ = li:,v^m +F ^- F ^ (16) 

, = Q f x dy fy 3 ^y 2 | g y\ 
Q Jo y — y^(y) \x 3 x 2 x J 



and hence the integrodifferential equation 



a r^{x-yfdy 

2 Jo 1 - EM (18) 



4 Aw \ --[ x i^yr^y 

X dx [X) ~2Jo l-E(y) 
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which likewise makes no reference to the subtraction point. We require the solution to 
that is regular for < x < fx 2 , with S(/x 2 ) = 0. Then the scalar anomalous dimension is 
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dlog(l-E(g 2 )) 



d\og(q 2 ) 



(1 — r) 2 dr 



2 Jo l-E(^r) 



(19) 



2.2 Multiloop integers from nonlinear differential equations 

By differentiating fll4] ) once and (|18|) three times we obtain 

(1 - t{q 2 ))P{D)t{q 2 ) = a = (1 - £(g 2 ))P(D)£(g 2 ) (20) 

with polynomials of the differential operator D := d/dlog(q 2 ) given by 

P{D) := D(D + 2) (21) 
P(D) := P(P + l)(P + 2)(P + 3) (22) 

and already encountered in analyses of rainbow and chain contributions. The rainbow 
diagrams generate contributions to the anomalous dimensions that solve 



P{ 7rainbow) — a — P\ 7rainbow) (23) 



7rainbow = 1 - Vl + « = ~~ + ^ ~ 2^ + 5^ + • • • (24) 



namely g, g 



3- \ 5 + AVl + a a a 2 a 3 

7rainbow = V - = -- + 1% - 206- + • • • 25 

2 6 o' 3 6° 

whereas insertion of chains in the one-loop diagram gives |l 

r°° , , xdx a a 2 a 3 a 4 



7chain 
7chain 



f°° . , xdx a a a a , . 

-2/ exp -2x/a ^ ~ — + — - 2— + 6— + • ■ ■ 26 

Jo V 7 ; P(x) 2 2 3 2 5 2 7 v ; 

f°° , . .xdx a a 2 „ a 3 , „. 



At n < 3 loops, rainbow and chain contributions correspond to same rooted trees. At 
n = 3 loops, there are 2 rooted trees and the full Hopf algebra gives the sum of the rainbow 
and chain contributions, with G% = 2 + 2 = 4 and G3 = 206 + 170 = 376. At n > 3 
loops, there are essentially new terms in the Hopf algebra. For example, G4 = 27 > 5 + 6 
receives contributions from all 4 rooted trees with 4 nodes. 



Thanks to the differential equations in ( p0[ ) we can achieve a dramatic speedup of 
the dimensionally regularized methods in |1[ 0. The ra-loop contribution to S((? 2 ) has the 
form a n J2k=i C( n > k)L k , with L := log(g 2 / /x 2 ). Substituting this Ansatz into the nonlinear 
fourth-order differential equation a = (1 — S(g 2 ))P(P)S(g 2 ), we require that all powers 
of the log disappear on the r.h.s. and that only the first power of the coupling survives. 
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Then the integer G n = 6(— 36) n_1 C(n, 1) is obtained by iterative solution of a system of 
equations that gives n linear combinations of {C(n,k) | 1 < k < n} in terms of 0(n 3 ) 
products of {C(m, k) | 1 < k < m < n}. 

This Dyson-Schwinger algorithm may be implemented by the REDUCE procedure 

procedure Hopf Scalar (n) ; for m:=l:n do 
<<s m:=for k:=l:m sum c k*log x~k; 

d m:=if m=l then 1 else for k:=l:m-l sum d k*s(m-k); 
sol:=first solve coef f (x*df (x~3*s m,x,4)-d m,log x) ; 
s m:=if m=l then log x/6 else (s m where sol); 
write G m: =6 (-36) ~ (m-1) *sub(x=l ,df (s m,x))»; 



which invokes the procedures solve, coeff and df , to solve for the coefficients of the 
results of differentiations. It yields the 30-loop coefficient (H) in 6 seconds, running version 
3.7 pi of reduce on a 500 MHz alpha machine. By way of comparison, the dimensionally 
regularized 30-loop calculation in U took 8 hours, on the same machine. The origin of 
the speedup is clear: here we have only a double expansion, in the coupling and the 
logarithm, whereas in [[l], § we also had a Laurent series in e, where the limit e — > was 
admissible only at the final stage of the dimensionally regularized BPHZ procedure. Now 
that we have a nonperturbative Dyson-Schwinger formulation, we can dispense with the 
BPHZ complexities. 

Note that, even using the new method, the time to reach n loops grows rather rapidly 
with n. At each iteration with m < n we solve a system of m linear equations and 
hence perform a number of elementary operations that grows like m 3 . The main cost 
comes from 0(m 3 ) multiplication of integers with 0(m log m) digits at the m-th step of 
the iteration. It follows that the time to reach n loops scales like n c , modulo logarithms 
of n, with an exponent c > 5, in any implementation, and with c < 6, in any sensible 
implementation. To achieve the lower bound c = 5 one might need to use fast Fourier 
transforms. We devised reduce code, more structured than that above, whose timing 
scaled, empirically, with an exponent c ~ 5.3, for n = 0(100). After a couple of days 
we reached 420 loops and checked the integer G i2 o by independent code that used David 



Bailey's mpfun package flTT |. As might have been expected, Bailey's multiple-precision 
extension of FORTRAN was somewhat faster. Table 1 gives the result for G500, obtained 
by running MPFUN for 21.5 hours. 

We now chart the discovery of a second dramatic speedup, which reduces the number 
of multiplications of large integers at the m-th step from 0(m 3 ) to merely 0(m). 

2.3 Parametric solution of the Yukawa problem 

For the Yukawa problem we work with the variable z := {q 2 /fJ?) 2 and define 



G{z) := ^ifa [z - (/iV /2 )) (28) 
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normalized so that G(l) = y2/a. Then the differential equation is simply 

2GG" = -1 (29) 
where primes denote differentiation w.r.t. z. This is easily integrated, to give 

(G 1 ) 2 = - log G + constant . (30) 
The next step is to work with the parameter p := G', in terms of which 



G = yj2/aexpW-p*) (31) 

where po is the value of p for which z = 1 and hence q 2 = /i 2 . Then a parametric solution 
is obtained by determining 

5(p) := z/G (32) 



which satisfies the first-order equation 



p 2p dp 

whose iterative solution gives the asymptotic series 



a = - H — (33) 



1 1 (2n-l) H 

p P ^ (-2p*y 



«(p) = : + :Etw- ( 34 ) 



From (|32|) and the derivative of fl28|) at z = 1, we determine 



a (Po) = \hi Po := -7= (35) 
V / V2a 



and hence obtain the asymptotic series @ at p = po. 

As advertised in the introduction, we can now give the complete nonperturbative 
solution of the d — 4 case, by noting that the asymptotic series (0) is that of 

roc 

5(p) = y/nexp(p 2 ) erfc(p) = 2 / <isexp(p 2 — s 2 ) . (36) 

Jp 

Tidying up, we then obtain the parametric solution 

~ 2 x \ja/2n 

exp(p 2 ) erfc(p) 

i = ^(^4) " 2 (38) 

with the condition (|l|) determining the anomalous dimension, even in the strong-coupling 
regime g > g cr i t '■= 2(27r) 3 / 2 « 31.5 where a > 2tc and hence : = (7 + 2)/\^2a < 0. 
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Finally we obtain the elegant recursion (^) for the coefficients of the perturbation se- 
ries (H), by turning (|33|), at p = po, into a nonlinear differential equation for the anomalous 
dimension. The result is conveniently written as 

27 = -d + a 2 |i- (39) 
da a 

which proves @ and generates the 500-loop Yukawa coefficient in 10 seconds. We remind 
the reader that this sums all diagrams at the 500-loop level generated by our Dyson- 
Schwinger equations, of which there are as many as there are rooted trees with 500 
vertices. Finally, we note that at g=30 the nonperturbative result for 7, to 120 digits, is 

-1.8520278058795936576120591680015771764936420535922291037938512550071476 
394977783846264445624911569322962723432940741627702, 

obtained by a Newton-Raphson method and in excellent agreement with resummation 



techniques 11 . 



2.4 Asymptotic parametric solution of the scalar problem 

In the scalar case, at d = 6, we work with the variable y := q 2 j \J? and define 

G(y):=j6/a~(y-yZ( t i 2 y)) (40) 
which transforms the fourth-order equation for E to 

gg " + U( 5+ 4) vG " = - 1 (41) 

where primes denote differentiation w.r.t. y. Taking p := G' as the parameter, we define 

a(p):=y/G; (3(p) := —GG" (42) 

with the formal solution 

generalizing (|36D , which was the Yukawa case (3 = 1/2. The differential equations 

1 (3 da . AA . 

a = - + --r G 44 
p p dp 

» = i+ il( 5 -0^ < 45) 

may be solved iteratively, to yield the asymptotic series 

aM ^ HSfw (46) 

«*> - >^,4 (47) 
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where A n and B n are positive integers, with A\ = 1 and Bi = 5, obtained by inspection, 
leading to A2 = 23 and B2 = 4 x 43, and so on. Diligent symbolic programming reduces 
the main computational burden to 0(n) multiplications of large integers, at the n-th 
iteration of 



By this means we developed 500 terms of the asymptotic series 



7-6£A n -— ^— (48 



n>0 



(67 + 6) 2 



in a minute, while a day was needed to compute G500, by the double-series method. 



2.5 Differential equation for the scalar anomalous dimension 

In the Yukawa case, we succeeded in computing G n via 0(n) large- integer multiplications 
in (|3|). Similarly, in the scalar case, we programmed the iterations of ( f4"4| , fS)D so as to 
obtain A n via 0(n) multiplications at the n-th step. 

Yet there is a glaring discrepancy between the Yukawa and scalar analyses, thus far. 
The procedure for obtaining G n from A n appears to involve 0(n 3 ) multiplications. Thus 
we have not yet improved on the double-series method, in the scalar case. We eventually 
remedied this problem, by a sequence of computer-algebra explorations and finally reduced 
it to the procedure DysonScalar to which we now turn. 



First, we sought generalizations of the method of ( f4^f45|) . The factors of (|22|) reveal 



that so far we have exploited only one of three possibilities. More generally, we obtained 

1 3 k da k 16^ A nk 
p p dp P P n> \—ujr) n 

1 + ^ip ( m ~ 6 _ eaA Tp) aA (50) 

(51) 




n>0 



(67 + 6k)' 



which is valid for k = 1, 2, 3, but for no other value. The analysis in (53,J3,^S|) corresponds 
to the case with k = 1. 



To check that (|50D carries the same information for k = 1,2,3, we used computer 
algebra to transform it to a nonlinear third-order differential equation for 7, invoking the 
constraint (k — 1) (k — 2) (k — 3) =0 to obtain 

8a 3 7 {7V + 477V + (Y) 3 } + 4a 2 7 {27(7 - 3) 7 " + (7 - 6)( 7 ') 2 } 
+ 2a 7 (2 7 2 + 67 + 11) 7 ' - 7(7 + 1)( 7 + 2) (7 + 3) = a (52) 

where the primes denote differentiation w.r.t. the coupling a. This seemingly baroque 
equation then remedies the previous discrepancy between the scalar and Yukawa cases, 
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allowing an iterative determination of the coefficients G n via 0(m) large-integer multipli- 
cations at each m < n. The iterative procedure is based on the rearrangement 

67 = -a + lla 2 ^-— -6T3 + T4 (53) 
da a 

T 3 := 4a 2 7 {77"+ ( 7 ') 2 } -2a 7 2 7 ' + 7 3 (54) 

T 4 := 8a 3 7{7 2 7 , " + 477V / + (7 , ) 3 }+4a 2 7 2 {277 ,/ + (7 , ) 2 } + 4a7 3 7 , -7 4 (55) 



where the first two terms on the r.h.s. of ( p3|) parallel those in (|39| ) and the remainder are 
grouped in T3, which is cubic in 7, and T4, which is quartic. 



At first sight, it was hard to spot the pattern of higher-order terms in (|54]j55D . It 
emerged when we discovered that 

T 3 = 7 (2aA-lW T 2 :=a^ 1 - = 1 (2a±- -lU. (56) 
\ da J dx a \ da J 

It was then natural to suppose that 

T 4 = 7 (<2aj- a - T 3 (57) 
and to verify that this indeed gives fl55l) . Hence we achieve a highly systematic iteration 



procedure DysonScalar (n) ; for m:=l:n do 

«for j : =1 : 3 do g( j + 1 ,m) : =f or k:=j:m-l sum g(m-k) * (2k-l) *g(j ,k) ; 
write G m:=llg(2,m)/6-6g(3,m)+6g(4,m)+if m=l then 1; g(l,m):=6G m»; 

which generates the 500-loop integer of Table 1 in merely 7 minutes. We will derive 
this procedure from the Hopf algebra structure of the perturbation series X, but first we 
exhibit a further consequence of this iteration. 



3 Momentum dependence via the anomalous dimen- 
sion 

3.1 Propagator— coupling duality 

Finally we arrive at a fine duality, between the dependence of the propagator 

nfe2) := (58) 

on the logarithm L := log(g 2 //i 2 ), and the dependence of the (suitably normalized) cou- 
pling constant y/a := g/(47r) d//4 on the dual of L, say A, which we shall now identify. 
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First consider the Dyson- Schwinger equation for II(g 2 ), which may be written 

4 2 ^)n^ = - n(AV5 - (59) 

Now consider the apparently baroque third-order equation (|52|) for the anomalous dimen- 
sion. Its origin is far simpler: 

p (~ 2a ^l) 7=a = (60) 

with the same polynomial, of a different operator. Thus the duality is as follows. 

Propagator-coupling duality At fixed coupling do let 

U(q 2 )q 2 ^T = S(L); L := log(g 2 /^ 2 ) • (61) 
Let 7(a) be the anomalous dimension at arbitrary coupling a. Then 

^- f a o db 

^ S(A > ; W. (62) 

Thus we expose the physical reasons for the enormous speedups that we have gained. 
The triple expansion in a, L and e, employed in |l], @, is obviously reducible to a dou- 
ble expansion in a and L, since we dispense with dimensional regularization, by taking 
derivatives of the Dyson-Schwinger equation. What was not so obvious, at the outset, is 
that the double expansion in a and L, underlying the procedure Hopf Scalar of sect. 2.2, 
is also unnecessary. Now, thanks to the amazingly compact procedure DysonScalar of 
sect. 2.5, we have a procedure for the perturbation series of 7(a) by iterative solution of 
the momentum-independent equation 



We had a strong intuition of the possibility of this second speedup, by a huge factor of 
0(n 2 ) at n loops, arguing as follows. Suppose that we knew the momentum dependence 
of the self energy, S(g 2 ). Then the anomalous dimension would be known, from the 
derivative of the self energy at q 2 = fi 2 , where S(^ 2 ) = 0. Conversely, one expects that 
a knowledge of the dependence of 7(a) on the coupling a is sufficient to determine the 
momentum dependence of the self energy, and hence that it is unnecessary to expand in 
both a and log(g 2 /yU 2 ). The dictionary (|6l| , |62|) establishes this fact. 



In summary: the self energy must be renormalized; the anomalous dimension encodes 
the infinities that have been subtracted. From the momentum dependence of the self 
energy one can find the anomalous dimension; in massless theories the converse is true. 
Note that we can read the duality as saying 

log(g 2 //, 2 ) f 

or as 

Theorem 1: 

dlog(l - 



d log q 2 

We now derive this result rigorously. 



q 2 n(g 2 ) x 7(x 2 a) 



7(0/(1 -s) 1 
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3.2 Proof of the duality 



We have to prove Theorem 1. We will rely on two propositions. The first exhibits how 
the Hopf and Lie algebra structure of perturbation theory interferes with the equation of 
motion -the Dyson-Schwinger equation- which is the defining equation X = aB + (l/ (1 — 
X)) for the Hopf algebra element X. The second proposition establishes the recursion 
DysonScalar. 

Let X = X^fcli ak Xk be a formal series in the Hopf algebra H of undecorated rooted 
trees , such that the counit e vanishes, e(X) = 0, with X k e H and a be the coupling 
constant, and as announced above let X be defined by 

X = a5 + (1/(1-X)), 

where we write 1/(1 — X) = X)£L X J '. This determines the X k e H uniquely. The first 
four X k read 

X 1 = B + (e) 

X 2 = B + {X X ) 

X 3 = B + (X 2 ) + B+(X 1 X 1 ) 

X 4 = B + (X 3 ) + 2B+(X 1 X 2 ) + B+(X 1 X 1 X 1 ), 

which are linear combinations of rooted trees, and the unit e of the algebra represents the 
empty tree. 

Note that the weights in X4 are not the Connes-Moscovici weights, but field-theoretic 
weights. 

Let Ft(U) be the befooting operator (with F Ul the usual befooting operator of |], |2"Tf ) 
which sums over all ways of removing the tree T from U, extended by linearity: for 
Y = cTi + dT 2 a linear combination of rooted trees (ie an arbitrary element of the linear 
basis of H), let us set 

F Y {U) = cF Tl {U) + dF T2 {U). 

Prop.l: The above coefficients Xj of the series X fulfill 

F Xm {X k ) = [2{k - m) - l)]X fe _ m , k > m > 1. 



Proof: The structure of the series X gives the following form to the coefficients X k : 



Xk+i — B + 



E E x%...xi 

0<i 1 <...<i r 



Si! . . . S r 



Vfc > 1. 



(63) 



Note that by definition Fx m {X k ) = 0, m > k and that Fx m is a derivation, Fx m {UV) 
F Xm (U)V + UF Xm (V). Also,™' 



F Xm (X m+l ) = Xi, Vm> 1, 



(64) 
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as 



Fx k {X, 



k+l> 



Fx k 



B A 



\ 



E E 

r = l siii + ... + s r i r — k 
0<il<...<ir 



Xf r - 



sil 



F Xk (B + (X k )) = X,. 



For any m in Fx m {Xj), we assume the proposition holds for all Xj, % < k. We show that 
it holds for X k+1 , Vm. For any m, Eq.(^) provides the start of the induction. 

Now, Fv and B + have the commutator 

[F Xm ,B + }(X k )=5 m , k B + (e), 

by the definition of Fx m {Fx m removes feet, B + grows at roots, only at m = k this 
interferes in the indicated manner), which generalizes to any product YliX ri of elements 
X Ti so that the interchange of the befooting operator generates some derivatives of this 
product: 



F x m B^ 



Y[x r 



B A 



^ m (n^)+E^n^ 



where the superscript fj at the product indicates an omission of the factor X T . . 

We can now interchange the action of the B + and the Fx m operator in the expression 
for Fx m (Xk+i), and use the assumption to find the desired result. Let us work out an 
example before we give some more details. 

X 6 = B + [X b + 2X : X 4 + 2X 2 X 3 + 3X^X3 + 3XxX 2 2 + 4X^X3 + X\}. 

For F Xm (X 6 ) we thus find (using F Xm {X x ) = 0,Vm) 

B + [F X JX 5 ) + 2X 1 F Xm (X 4 ) + 2F Xm (X 2 )X 3 + 2X 2 F X JX 3 ) + 3X 1 2 F Xm (X 3 ) 
+GX 1 X 2 F Xm (X 2 ) +4X 1 3 F Xm (X 2 )] + commutator terms. 

To be specific, let us set m = 2 to get, using the assumption of the induction, 

5+ [5X 3 + 6X!X 2 + 2X 2 Xi + 3X 2 Xi], 

while the commutator terms deliver 

5+[2X 3 + 6X 1 X 2 + 4X 1 3 ]. 

Altogether, we get 

£? + [7X3 + 14XxX 2 + TXj 5 ] = 7B+ [x 3 + 2X X X 2 + Xf] = 7X 4 , 

as desired. 

Let us now turn to the general case. Let us consider a monomial of degree r contribut- 



ing in the sum in Eq. 



It suffices to consider the generic case si 



1, if 
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there are higher powers Sj > 1, the fact that Fx m is a derivative ensures that the same 
argument holds and that the extra factors in the denominator will be appropriately 
cancelled. 

In the generic case the coefficient of the monomial is r\. Let us first consider the case 
r = 2. Contributions to X a X b say, with a + b = k — m, are coming from 

2*Fx m (B+(X a+rn Xb + X a X b+m )). 

By assumption of the induction, this delivers, for the non-derivative terms, 

2\[2(k-m)-2]B + (X a X b ), 

while from the derivative terms, we get 

3\B + (X a X b ), 

which adds up to 

2\[2({k + l)-m)-l)B + {X a X b ), 

as desired. In general, for arbitrary r, the simple fact that r!(r + 1) = (r + 1)! en- 
sures that each term in the sum homogenously factors [2((k + 1) — m) — 1], and sum- 
ming over all monomials and applying the B + operator, we get, by Eq.(|63|) backwards, 
F Xm (X k+1 ) = [2((k + 1) - m) - l]X fc+1 _ m . □ 

Now, let <7^ := lim e _, e n [0 o [Y 11 ]] (X m ), with antipode S and grading Y and convo- 
lution product * as in [[|. 
Then, the procedure DysonScalar is 
Prop. 2: 

m— 1 

£[2(™-i)-iW<ci- 

8=1 

Here, <fi is the unrenormalized character on the Hopf algebra <fi = X |„2 =1 . With its help 
and the Hopf algebra automorphism Q p : H — > H of ||19|| , 



6 p (/i) = e pdcg{h) h, VpeH 

(deg is the degree: deg(X m ) = m) we can write the renormalized character: the character 
S is given by 



-elog(<j 2 /^t 2 ) 



The Taylor expansion in log(g 2 //i 2 ) is then solely a study of the expansion of the Hopf 
algebra automorphism S * 6_ e i og ( g 2/ /J 2) which can be written, setting z = log(g 2 //i 2 ), as 

S * 0_ e i O g(g2 /At2) = — — h~ ez Y- 
i=i % - 

Also, we note that quite generally 

e n (f>(S * r[I 3 ]) = C(e n " m ), VXj, (65) 
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which follows from the fact that S exists in the limit e — > 0. 

Proof: We use Prop.l and induction on m. We only consider the case n = 2, the other 
cases follow identically. The start of the induction involves an explicit check of properties 
of 4>(Xi), (f)(X 2 ), which can be done easily using 4>(Xi) = 5 1 ,0(X 2 ) = BiB 2 , where 

Bl = I dDk [k^(k + g y ' 92=1 ' 

We concentrate on the induction step. Assume 

e 2 <j>[S * Y 2 [X 3 }} = e 2 £ <f>[S * YfaMS * Y[F Xi (X,)]] 
i=i 

for all j < m. We want to prove it for m + 1. Again, we use the presentation Eq. (|63|) 
for X m+ i, for which we write in shorthand X m+1 = B + (X m + J2Hi^i)i empasizing a 
decomposition into product terms and the single term X m which is linear in generators 
in the sum involved in this equation. 

We then have 

e 2 <J)(S*Y 2 [X m+1 }) = e 2 (j)[S*Y 2 [B + {X m + Y.Tl X i)\\ 

i 

= ^s^l^WD + ^Ms^l^Il 1 !' ( 66 ) 

i 

Here, B + (</>) = <p o B + and we use 

S*Y 2 oB + = mo(S®Y 2 )o(id® B + ) o A 

(as Y 2 (e) = and where m is the product, A the coproduct) which further can be written 

as 

m o (id ® B+) o(S® Y 2 ) o A, 

due to the fact that for the commutator between Y,B + we have [Y, B + ] = B + . Hence, 
interchanging B + with Y 2 produces only commutator terms involving e 2 (f) o S * Y l , with 
% = 0, 1, which vanish by Eq.(^) in the limit e — >• . 

Let us treat the two terms on the rhs in the last line of Eq.(^) separately. For the 
second term, we use that S * Y k [l\ i X- h ] = if the product has more than k factors. To 
see this, consider, for example (using Sweedler's notation), 

Sieves) = SiXlX^Xi'X^degiX^+degiX'!)] 

= S * id(Xi)S * Y{Xj) +S* Y(Xi)S * id(Xj) = 0, (67) 

and similarly for higher powers of Y. Hence, S * Y 2 vanishes on products of more than 
two arguments, and on two arguments we have, S * Y 2 (XiXj) = 2S * Y(Xi)S * Y(Xj). 
(In general, S * Y n behaves like a n-th derivative acting on a product of functions which 
vanish at the origin, for example 

(fg)" = f"g + V'g' + fg" = V'g'-) 
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Hence, the second term on the rhs above delivers 

2e 2 £ B^ttSKYlXiMSKYlXj]]. 

i+j=m 

For the first term on the rhs we find 

e 2 [<p ® <f)](S ® F 2 )(id ® 5+)A(X m ) = m[(id ® S^)(0 ® <8> F 2 )A(X m )]. 
Here we interchanged B + with the K operator as before. The assumption of the induction 



is 



{S*Y 2 ) o A(Xj 



S *Y ® (f>o S *Y) o A(Xj). 



From this assumption and from Prop.l we can write for the above expression 
Y\2{m - i) - 1] [0 o S* Y{Xi)\ [B + {4>) o S * Y(X m -„ 

i 

Adding the results of the first and second terms for the rhs of d66|) , we find 

£[2((m + 1) - i) - 1] [0 o S * KpQ)] o 5 * F (X m _,) 

j 

Using that S"*y vanishes on products by Eq.([)7]), we can combine the resulting expressions 
to the desired result, by adding zero in form of adding all the product terms in Eq. flBlf ) 
for X m+ i_j to eliminate the B + operator in favour of the increased degree m+l — i. This 
proves the proposition. □ 



With these tools, the proof of Theorem 1 becomes a mere comparison of coefficients of 
two series in an infinite number of variables. Indeed, in both series we can express any 

m. Hence, Theorem 1 



re- 



appearance of cr^j by Prop. 2 as a product n™ =1 o"^ , 
provides two power series in an infinite set of variables a\. Note that a monomial in 
degree m in these variables is always accompanied by the factor z" 1 " 1 . 

Let us set 



oo oo n „n 



EE 

k=l n=l 



111 



and 



a 



E 

k=l 



Then, Theorem 1 amounts to showing that 



oo oo r 



1 



k=l n=0 



nl 



1 



a 



oo 

E a m 

m=l 



(I -a) 



2m 



(68) 



Using Prop. 2 and elementary algebra one finds the following coefficients at orders 
z°, z 1 , z 2 /2, ... for the lhs: 



a 



(i) 
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z 1 : a^^i^K, 

ri=l 

7 2 oo 

-: Y {2r 2 ){2r 1 + 2r 2 )al i a^ 

n,r2=l 

which for example implies that the coefficient of the term [u\\ n z n ~ x / (n — 1)! is simply 
2 n ~ 1 (n — 1)!. The crucial step in all this is to recognize how the terms coming from the 
expansion of 1/(1 — a) cancel the ' — l'-term coming from Prop. 2, ie the term — J2i cr l cr '^i 
in 

< = E[2(m - i) - iK x <r-V 

i 

A similar phenomenon happens in the Taylor expansion of 

1 

(i-Eti^ 1 ) 2 " 1 

to the fc-th order, involved in the rhs of (pip. For example, at order z 2 , k = 2, the ' — l'- 
term in (2m — 1) generated by the second derivative does the same job, and in general 
the lhs and the rhs of ( |S8"D give the same series in z. □ 



4 Conclusions 

Returning to the questions in the introduction, we answer them as follows: 

1. The non-perturbative problems whose expansions were studied in [|IJ Q are defined 
by (14,18). For the Yukawa problem, a solution in closed form is given by (1), for the 
anomalous dimension. The duality of Theorem 1 then gives the propagator as a limit of 
integration of an integral that must evaluate to log(g 2 //i 2 ). The particular solution of this 
duality is given by (37,38), in terms of the complementary error function. The same form 
of duality applies in the scalar case, though for this we have not yet obtained an explicit 
solution to the 4th-order equation that results from (53-55). We have, however, found 
a way of developing the perturbation series almost as efficiently as in the Yukawa case, 
attaining 500 loops in a few minutes. 

2. The obstacle in (5) at 7 = —2, i.e. at coupling g = 2(2n) 3 ^ 2 rs 31.5 in the Yukawa 
case, was surmounted, thanks to the parametric non-perturbative solution (1). In the 
scalar case, the obstacle in (48) at 7 = —1 is surmountable, by numerical solution of the 
third-order equation for the anomalous dimension. In the Yukawa case, we are able to 
compute strong-coupling results to very high accuracy and use them as stringent tests 
of the Borel resummation in |TTJ; in the scalar case, we have confidence in such strong- 



coupling resummations, but cannot yet surpass them with the ease afforded by (1,37,38). 

Let us now finish this paper by a discussion of the difficulties one confronts in general. 
Our starting Dyson-Schwinger equations were addressing Feynman graphs representing 
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only the Hopf algebra of undecorated rooted trees. To treat the general case, the first 
observation is that there would be various operators B + , indexed by labels corresponding 
to primitive graphs without sub divergences, of the form 

W * _ a deg{x) r dH F*(l 2 J-q,q 2 ) 



7r 2deg(x) J 12 ' 

where deg gives the loop number of the decoration x, F x is some integral kernel obtained 
from the skeleton graph x, and the corresponding integral operator contains a sum B + = 
J2 X B+ ■ Here, the previous representation of the B + operator as 

a f dH _ q-l 



it 2 J l 2 (l + q) 2 

would be just the first term in such a series. 

A further complication is given by the fact that allowing for internal self-energies 
and vertex-corrections one would have to consider a coupled system of Dyson-Schwinger 
equations for those functions. Still, all this is conceptually not too drastic a challenge, 
and Hopf and Lie algebra structures for the full theory are readily available || [18[ . 



But the real difficulty resides in the fact that these new B + operators offer a choice 
of as many places as there are internal lines and vertices for the self-energies and vertex 
corrections to be inserted into x. This is in sharp contrast to the Dyson-Schwinger equa- 
tion studied here which was carefully designed so that it was always the same propagator 
into which the self-energy was inserted (see Figure 1). 

It would be very convenient if the characters of the full Hopf algebra of Feynman graphs 
would turn out to be independent of such choices, alas, in dimensional regularization, they 
are not. Nevertheless, they fail to be so in a most interesting way, as mentioned in f20|| : 



only lower order pole terms are affected by such choices, and a proper understanding of 



the operad structure of Feynman graphs |20| should ultimately allow to generalize the 
methods presented here to the full theory. 
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Table 1 Coefficient of the scalar perturbation series with 1675 digits at 500 loops 

G500=206261451966080541451119356265266407905816117576895601520616328670543304097 
62369668214104674763068056454522518617422020409397336434904863988900797769773644 
47129884863324773181376863120291798830884688213932683869821267125662274428136514 
68974978228592824043044373847281757207937081063432528806815509319762088807291996 
54549245884853496719417048678199825379018355919198123075612308008976364608893906 
00835837012056033720017238115336850340799075684336975651857656078799282745256216 
85768456030809283727097722850488278232311177219444745322287340871435443707536590 
64304859950724683157717734493071321199539578218428617617722892100276682781401203 
04983974209704793621909710059353724523231635766062166284812903992269403282699432 
81718327508638643305481989940132234093616573076862094588977827344981584305605437 
66475002382217933275761312682929603923397580260987048907414858143897114762331252 
08694985337972553885925402003826420205441859988844001867088083850782378303677991 
14077650584544145709672328391394562704209221732180879565868213522109303655045186 
92714017665002971967455255310508358729281544729403249398746232320441525286283859 
23093041626365262630100048817481274793707664791767175677240144896307853488347045 
21622394885797995125083750860330519417878429051836575477220881369445751634601965 
33191009573619480068718080810533581305996863996579338522874547127421808710757882 
86996199556804886954946559116947132125235605586627322129268965041445488085748194 
82341875039156647569797757032552836429751077302524927736861138479038542006096835 
73747720303607608007740173613335602076396299832459826245418033598839559699294537 
37336134624690115674194793212055897162647586497730033948880084738561472545509216 
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